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Abstract 

Q_j The many-body wave function of a pion and its decay products shows that the probability 

of detecting a neutrino possesses a unique finite size correction. The detection rate at a finite 
distance L is expressed as Tq + g{ujuh/ c)Ti, where oj^ = m'lc'^ /2Eyh and c is the light velocity. Fq 
is a constant that is computed with the standard S-matrix of plane waves and the second term is 



04 a finite size correction that is computed with wave packets. The value of g{x) decreases rapidly 

7^ ; with X and vanishes in charged leptons, but is finite in neutrinos at a macroscopic L. The finite 

size correction reveals a diffraction pattern of a single neutrino and can be computed rigorously 
with the light-cone singularity of a system consisting of a pion and a muon. We predict that the 
neutrino diffraction would be observed at near-detector regions of ground experiments and that it 
could be used for the experimental determination of the neutrino mass. 
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1 Neutrino interference. Interference phenomena of photons, electrons, neutrons and 
other heavy elements are useful for testing quantum mechanics and other basic principles. 
In this paper, we show that a neutrino produced by pion decay displays an interference 
phenomenon that could provide an absolute value for the neutrino mass from its unique 
interference pattern. 

A neutrino interacts extremely weakly with matters; hence, a quantum-mechanical wave 
function that expresses a system consisting of a pion and a charged lepton, and a neutrino 
wave, become highly correlated. Consequently, the detection probability of the neutrino 
becomes dependent on a distance between the initial pion and the neutrino detector, which 
we call a finite size correction. An S-matrix ^[T] defined according to the boundary condition 
at the finite size, where T stands for a time interval, does not satisfy various relations [l| of 
an S-matrix of the infinite size S'[oo]. Wave packets, localized around center positions, satisfy 
the asymptotic boundary conditions I2, 3| even in this case and are suitable for studying 
the finite size corrections of transition probabilities. Accordingly we compute the transition 
amplitude with the S-matrix of the finite time interval, 5'[T], defined using wave packets. 

We find that the detection rate is composed of a constant and a finite size correction. 
The former agrees with a standard value obtained by an S-matrix of plane waves, while 
the latter is a new term that cannot be computed by it but by ^[T] and has an origin 
in a diffraction. Waves accumulating at the light velocity form a light-cone singularity and 
exhibit a neutrino diffraction phenomenon, which makes the finite size correction long range. 
The diffraction pattern is determined by the difference of angular velocities, = — w^^, 
where = Ey/h and ufP = cIpMh. ujy becomes an extremely small value mlc^ /2Eyh 
owing to unique neutrino features |#^. Consequently, the diffraction term becomes finite 
in a macroscopic spatial region r < "^^J^^Jr and affects experiments in a mass- dependent 
manner at near-detector regions. 

A neutrino wave packet jr-Q] expresses a nuc 



the neutrino interacts and is well localized 



ble {4], thus, we study a situation in which the mass-squared average ml satisfies, ^ ^^l-, 
and present one fiavor case. Extensions to general cases are straightforward. 

2 Position-dependent probability. The position-dependent probability of detecting 
a neutrino is computed with S\T], which has two components owing to the algebra of [1| 
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eon wave function in a nucleus with which 



16l ] . Mass-squared differences 6m1 are negligi- 



and is expressed as 

(/3|5[T]|a) = (/3|5(")[T]|«) + (/3|5W[T]|a), (1) 

where the first term satisfies E13 = Ea and the second term satisfies Bp 7^ Ea- The pro- 
cess by which a pion is prepared at one time and a neutrino is detected at another is 
described by the time- dependent Schrodinger equation and the position-dependent amph- 
tude T = f d'^x {fi, h'\Hw{x)\TT) between a pion prepared at T,r, a wave packet of a neutrino 
at (Tj,,Xj,) and an unmeasured-muon. They are expressed as |7r) = 1/4, T^), If',^) = 
I'yPu, Xy, T,^). T is written with the pion field and Dirac spinors as 

d^xdkjy N{0\if^{x)\7T)u{pf,){l - 

y ^ip^i-x+ikv-{x-Xv)-^(k,,-pv)^ 

where N = igm^ {a,,/7i)^ {m^nii,/ E^E^,)^ , and t is integrated in the region T^, < t. The 
value for must be integrated in the whole region in order to represent the localized wave 
packet in coordinate variables. For simplicity, a Gaussian form is assumed in most places. 
We will show that the finite size correction has a universal property that is common to 
general wave packets. The size of the wave packet, cr^, is estimated later. The amplitude T 
depends upon the coordinates (T,^, X,^) and T^ and the probability depends on T = T,^ — T^^. 

Integrating k^, we obtain a Gaussian function of x, which vanishes at large \x\ and 
satisfies the asymptotic boundary condition. We express the square of the modulus of T 
with a correlation function. Because the order of integrations is interchangeable, the muon 
momentum is integrated first for a fixed Xj. Then, after the spin summations are made, the 
transition probability is written with the correlation function A[xi,X2), 

= / ^fe.pje-'--'- (2) 

as, 

Eu J 
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where C = g^m'j^ {Air/a,^)^ V is the normalization volume for the initial pion, x° = 
'^u + Vuiti — T^), 6x = xi —X2, 4>{Sx) = pi^-6x. In Eq. (|2]), the muon momentum is integrated 
in the whole positive energy region in order for Eq. ([3]) to agree with the original probability. 



3 Light-cone singularity. ATr^fj,{Sx) has two components. One is singular and extended 
in large \6x\ and the other is regular and localized in small \6x\. A superposition of infinite 
waves of the same phase leads to the first component. This term is easily extracted by an 
expression of the integral in the four-dimensional form with the new variable q = ~ Pn 
that is conjugate to 6x. Then A^_^((5x) is decomposed into the integrals of the regions < g° 
and — < 0'° < 0. The integral in the region < g° is expressed as, {pn-Pu — 
where 



4 0{q'' 



Im 



+ 2pT,-q + — ie 



and m 



m: 



m^. The integrand of Ji is expanded in p.„-q and the integration leads to 



the light-cone singularity 



17 1 , 6{6x'^), and less singular and regular terms that are described 



with Bessel functions. The integral from the region — < 5''^ < 0, /2, is written with the 
momentum q = q + Pn and has no singularity. Thus the correlation function, At^^^{6x), is 
expressed as 

d 



d6x 



d 



+ f short + /: 



(4) 



d5x J \ Atx 

where A = {5x)\ D^{-i^) = (1//!) (2p.-(-^ai) ^)', fskort = -fj ^(-A) {N^iO' 
i e{5t) Ji{^)}— 9{\) Ki{^), ^ = rh\^, Ni, Ji, and Ki are Bessel functions, fshort has 
a singularity of the form 1/A around A = and decrease as e~™V^' or oscillates as e*"^V^I 
at large |A|. The condition for the convergence of the series will be studied later. 

Integration of spatial coordinates. Next, Eq.dl]) is substituted in Eq.(l3]) and xi and 
X2 are integrated. The first term, Js{\), is derived from the most singular term, ^^5(A), 
and computed with the center coordinate = {xi+X2)/2 and the relative coordinate 
r = xi — X2- Integrating these variables, we have Js{\) as 

e{6t) c^,(st)--^5t^ 



J5{\) = C, 



5(A)- 



(5) 



where C5(a) = (o'i/7r)2(j^/2 and (t)c{,5t) = Ui,6t = 5t mlc /2E^. We note that the phase 
(f){6x) of Eq. ([3]) became the small phase (pd^t) of Eq. at the light cone A = 0. The next 



singular term is from 1/A in Att^^^, and becomes Jst^X)/ ^T^a^pl, which is much smaller than 
Jsi^X) in the present parameter region. The magnitude is inversely proportional to \5t\ and is 



independent of m^. This behavior is satisfied in general forms of the wave packets and the 
following theorem is proved. 

Theorem: wave packets of spreading non-Gaussian forms In general cases, the 
leading term Js[x) of the wave packets of spreading non-Gaussian forms has the phase and 
the magnitude of the form Eq. The phase is given as a sum of (pdSt) and a small 
correction that becomes 0{1/Ey) in general and 0{1/El) in time inversion invariant wave 
packets. 

(Proof.) The Gaussian function of Eq.(l3]) is replaced with the wave function ip{x — vt) 

X ipkik, kr), Cij = 6ij/2E, (6) 

where the last term in the exponent is from an expansion of E{p + k) and makes the wave 
packet spread with time and ipk{ki, kx) is a wave function in the momentum. Since the 
coefficient Cij in the longitudinal direction is negligible for the neutrino, it is neglected. The 
center coordinate X is easily integrated in Eq. ([2]), and the product of ip{x — vt) and its 
complex conjugate is reduced to ipiri — v^dt^rx) and is written as 

i){ri-vJt,rT) = j dkidkre''"^'^-"-^''^ 

xe"'^-'''^+'^''T/^^^^'\^k{ki, ^t)!'. (7) 

Then the leading term Js(x) = / dre^'^'^^^'^'ipif — Vy5t)^^^6{\), becomes Js(\) = 7re^5tyg(l + 
7)^e^'^=(^*)(i+'5) at a large \5t\ region, where Wq = j 0)p, 5 = di/E + d2/2E^, 

7 = di/2E + (^2/2!) X (1/2E)2 - (1 - v,)^6t\ di = wq-^ J dkiki\Mki,0)\^ and d^ = 
Wo~^ J dkikf\iljk{ki, 0)1"^. In a wave packet of time reversal invariance, |w(/c;,0)p is an even 
function of ki, hence di = 0. The correction terms become negligible for E^ > 2 — 300 
[MeV]. ( Q.E.D.) Thus the theorem is proved. This theorem demonstrates that the neu- 
trino fiuxes for any collision are modified with the diffraction component. 

Next we evaluate the integrals of the regular terms of A^^^^. These terms are oscillating 
or decreasing rapidly with A and those oi r {] contribute to this. Hence, the spreading 
effect is negligible, and the Gaussian approximation is good. 

The first term, Li, is from f short in Eq.(|l]). In the space-like region A < 0, the asymptotic 
expressions of the Bessel functions at large \5t\ give Li = Ci\5t\~i e'^^" "P"'''' )st-a,,pl+irhy/\2v,.a,yp, 
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where Ci = (^^) ^ {4:VuauPu)~* ■ On the other hand, in the time-hke region A > 0, Li 
decreases with time as e~™^^V^. The second term, L2, is from J2, which is approximately 
the integral of e~*(^''~'^''~V^^^+™l")''* in g in a range Thus L2 is a steep decreasing 

function of 6t. 

4 Time-dependent probability. Finally we integrate ti,t2 in. 



J \0t\ Oy Tl \ Tl / 



(8) 



where Ni = ig'^rn?^Ti^a^{^pT,-py/Ey)V^^. In the above equation, and most other places, the 
low neutrino mass is neglected compared to fh?, Pn-Pu and ai,, except the slow phase (pd^t). 
The first term in Eq. (|8]) oscillates slowly with time 6t and the remaining terms oscillate or 
decrease rapidly. They are clearly separated. 

Here we study the series ^n(~'^P'^'P'')^7J. (a^)'^-^i' Eq- dH])- This series converges 
when the most diverging term, = E„(-2p,r-Pi/)"^ {-d^T ("^^)^ = {^10^T T^-'^i^)^ 
becomes finite. This is clearly seen in 2p.,^-pi, < rh? . At 2pT^-py = m^. Si becomes finite, 
and the value is expressed with the zeta function, (5/4) (m)^. Hence, in the 2p^-pj, < 
region, the series converges. Then the power series rapidly oscillates with a/I^^I as 5*2 = 
^tm^y\2v^avp„5t\(^i-''r^^ )^ Therefore the present method is valid in the region 2pT^-py < rh?. 
Outside this region, the power series diverges, and A^^^(5x) has no light-cone singularity. 
Then A^_^((5x) has only the short-range term. 

The integrand of the first term in Eq. ([H]) is oscillating extremely slowly with 5t and its 
integral of ti and ^2 in the finite T 



r dtidh-^e^^^'' = T(^KT) - tt), (9) 

Jo M 



slowly approaches constant, while g{ui,yT) satisfies ^^(c<;j,T)|t=o = —^u and ^(00) = 0. 
The last term in Eq.(|9]) is cancelled by the short-range term Li in Eq. ([8]). Here g{u,^T) is 
generated by the superposed waves of the form of light-cone singularity and its effect remains 
in a macroscopic distance of the order ^^Ifr- We call this the diffraction term. 

The last term in Eq. is f a/v / dtidt2L2{St) = TGq, where the constant Gq is com- 
puted numerically. Owing to the rapid oscillation in St, this integral receives contributions 
from the microscopic \6t\ region, and consequently Go is constant in T. This term does not 
depend on a^, and it agrees with the normal probability obtained with the standard method 
of using plane waves. 
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FIG. 1. Total detection rate at a finite distance L is given. The constant shows the normal term, 
and the diffraction term is written on top of the normal term. The horizontal axis represents the 
distance in [m], and the normal term is normalized to 0.8. The excess is seen in the distance below 
1200m. The neutrino mass, pion energy, and neutrino energy are 1.0 [eV/c^], 4 [GeV], and 700(A) 
and 800(0) [MeV], respectively. 

Integrating the neutrino's coordinate X^, in Eq. (|8]), we obtain the momentum- dependent 
probability. The total volume emerges and is cancelled with V^'^ of the normalization of the 
initial pion state. The total probability becomes 



where N2 = STgf^m^cTjy and L = cT is the length of decay region. The first term corresponds 
to S''''^''[T], while the second term corresponds to S''-"''[T]. At T — t- 00, the diffraction term 
vanishes and the probability P agrees with the value of the standard calculation using plane 
waves. At finite T, the probability has the diffraction component, which is stable under the 
variation of the pion's momentum. 

Next, we evaluate each term of Eq. fllOj) . In the normal term Gq, the energy and momen- 
tum are approximately conserved, and the product of the momentums is expressed with the 
masses p^-Pu = 'fh'^/'2. Integrating the neutrino's angle, we find this term independent of 
CTjy, which is consistent with the condition for the stationary state 1^. The transition rate 
agrees with the value of the ordinary method. However, the diffraction component, ^(T, u,^), 
is present in the kinematical region, |pi^|(ii^^ — \pn\) < Pn-Pu < from the convergence 

condition and is integrated in this region. This is slightly different from p^- Pu = m?/'^, 
hence it is impossible to experimentally distinguish the latter from the former region. We 
add both terms. The total probability thus obtained is presented in Fig.[T] for the mass of 
neutrino, rriy = 1 [eV/c^], a pion energy i?^ = 4 [GeV], and a neutrino energy E^, = 700 and 




(10) 
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800 [MeV]. For the wave packet size, in which the size of the nucleus having a mass number 
A, ay = A^/m1 is used. The value becomes = QA/m'^ for the ^^O nucleus and this is 
used for the evaluation. From this figure, we see that there is an excess of flux at the short 
distance region L < 600 [m] and the maximal excess is approximately 0.2 of the normal 
term at L = 0. The slope at the origin L = is determined by Uy. The diffraction term is 
slowly varying with both the distance and energy. The typical length Lq of this behavior 
is Lq [m] = 2Eyfic/m1c^ = 400 x Ey[GeV]/ml[eV'^ /c"^]. The neutrino's energy is measured 
with uncertainty AE^ in the experiments, which is of the order 0.1 x Ey. This uncertainty is 
100 [MeV] for the energy 1 [GeV] and the diffraction components of both energies are almost 
equivalent to those in Fig.[Tl For a larger value of energy uncertainty, the computation is 
easily made using Eq. ( |T0|) . Hence the diffraction component is observable if the absolute 
value of the mass is around 1 [eV/c^] using the near detector, but it becomes difficult to 
observe if the mass is less than 0.1 [eV/c^] using the muon neutrino. In the latter case, an 
electron neutrino may be used. 

The process described with ^[T] has the total probability Eq. ( ITOl) . In the same experi- 
ment, the detection rate of the muon, after neutrinos are integrated, has the same excess. 
Ordinary experiments of observing the muon, however, does not observe the neutrino and 
is described by another S'[T'], which satisfies the boundary condition for the muon and 
T' = T^ — Ttt is a time interval for the muon observation. The muon detection probability is 
computed with a free neutrino, and a probability, then, is expressed in the form of Eq.l JTOll 
with LOy ^ = rn?^c^ /2E^h. Since the muon is heavy, w^T' becomes very large and g{uj^1') 
at a macroscopic T' vanishes. Thus the detection probability of the muon is not modified, 
and it agrees with the normal term. The light-cone singularity is formed in both cases, but 
the diffraction is large in the neutrino and small in the charged lepton. 

The detection probability of the muon depends on the boundary condition of the neutrino. 
When the neutrino is detected at T, the muon spectrum includes the diffraction component, 
but when the neutrino is not detected, the muon spectrum does not include it. The latter 
is the standard one, and the former is non-standard, but may be verified experimentally. 

Here we compare the neutrino diffraction with the diffraction of light passing through 
a hole. In the neutrino, the diffraction pattern is formed in a direction parallel to the 
momentum along with the phase difference ujy5t of the non-stationary wave. Its size is 
determined by Uy, which is extremely small and stable with variations in parameters. Hence, 
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without fine tuning of tlie initial energy, tlie observation is easily made. In the light, the 
diffraction is formed in a direction perpendicular to the momentum with the phase difference 
(jj'^^St of the stationary wave. Its shape is determined by oj'^^ , which is large and varies with 
a change in the parameters. Thus the fine tuning of the initial energy is necessary for the 
observation. 

5 Summary and implications. The neutrino detection probability was computed 
with the S-matrix of the finite time interval S\T] and received a large finite size correction 
in addition to the normal term. The latter agrees with the known value obtained with 
the standard S-matrix of the plane waves, whereas the former varies with T in the form 
Cig{ujJT), where the constant Ci depends on the wave packet size. Here giuJT) slowly varies 
with T and the size of the diffraction pattern is determined by the small angular velocity 
Up instead of the de Broglie wave length and it is stable under the change in parameters 
of the initial states. Thus, a neutrino inside this distance reveals the quantum-mechanical 
diffraction of wave natures and this is not expressed by the asymptotic states of plane waves. 

The diffraction gives new corrections to neutrino fluxes but not to those of charged 
leptons, thus, it is consistent with all previous experiments of charged leptons. The neutrino 
diffraction in pion decays is a phenomenon distinctive to neutrinos in the non-asymptotic 
region. 

In the case of three masses mj,. and a mixing matrix Ui^ai the diffraction term to a 
neutrino of flavor a is expressed as 'Yl,i9i.'^i^i^)Wi,a\^ ^ whereas the normal term is expressed 
as I Ui^^D{i)U} where i is the mass eigenstate, a is the flavor eigenstate, and D{i) is 
the free wave of 711^^- Hence the diffraction term depends on the average mass-squared m^, 
but the normal term depends on mass-squared differences 5m^. At L — )■ 00, the diffraction 
term disappears and the normal terms remain. 

The new term has various implications for existing neutrino anomalies and future ex- 
periments. One anomaly is an excess of the neutrino flux at near detectors of ground 
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experiments. Fluxes measured by the near detectors of K2K [18'] and MiniBooNE 
excesses of 10 — 20 percent of the Monte Carlo estimations, whereas the excess is not clear in 
MINOS {20!. These excesses may be connected with the diffraction component. With more 
statistics, quantitative analysis might become possible to test the diffraction term. Another 



is the LSND anomaly 2l|] in which electron neutrinos in pion decays had excesses. Since the 



diffraction has the origin in the energy-nonconserving term, the helicity suppression does not 



work. An electron mode retains in the neutrino flux. Hence the excess in the near-detector 
regions is attributed to the diffraction component, instead of the flavor oscillation, and the 
controversy between LSND with others is resolved. 22| Finally, the distance or energy 
dependence of the neutrino flux may provide a new method for determining the absolute 
neutrino mass. 

Thus the neutrino diffraction appears visible at macroscopic distances and would be con- 
firmed with near detectors. At much larger distances than the above length, the diffraction 
component disappears and only the normal component including the neutrino flavor oscil- 
lation, remains. If the masses do not satisfy ml ^ 6ml but satisfy ml ^ Sml, then the 
neutrino fluxes show more complicated behaviors. 

A new quantum phenomenon of neutrinos specific to extremely small masses was derived, 
and its physical quantity determined by the absolute neutrino mass was presented. 

In this paper, we ignored the higher-order effects such as the length of pion life time and 
the pion mean-free-path in studying the quantum effects. We will study these problems and 
other large-scale physical phenomena of low-energy neutrinos in subsequent papers. 
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